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We summarize recent work on a frustrated periodic long-range Josephson array in a parameter 
regime where its dynamical behavior is identical to that of the p = 4 disordered spherical modeL 
We also discuss the physical requirements imposed by the theory on the experimental realization of 
this superconducting network. 
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The identification of the key features underlying the physics of glass formation remains an outstanding problem of 
statistical mechanics. A glassy system has a "memory" of its past behavior and does not explore all of phase space; 
it breaks ergodicity without thermodynamic selection of a unique state thereby defying description by the standard 
Gibbs methods. Though glass formation in the absence of intrinsic disorder is a widespread phenemenon, it remains 
poorly understood even at the mean-field level. Because vitrification is a dynamical transition that is not necessarily 
accompanied by a static one, it remains outside the framework of conventional Landau theory. The crucial links 
between glasses with intrinsic and self-generated disorder remain an area of active discussion. Recently several 
microscropic non-random models for glassiness have been proposed; most were studied via mapping to disordered 
systems. Q What common/distinct features are characteristic of glasses with and without randomness? How sim- 
ilar/different are glasses with short- vs. long-range interactions? What aspects of existing theoretical models are 
relevant for experiment? In this Proceeding we summarize our current efforts to address these questions in a periodic 
model that is both analytically accessible pM and experimentally realizable. |H| 
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Fig. 1 The phase diagram of the array (insert) where Tq indicates the temperature associated with the dynamical instability discussed 
in the text, Ts is the speculated equilibrium transition temperature and Tm is the "superheating" temperature where the low-temperature 
metastable states cease to exist. 
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The physical system under study is a stack of two mutually perpendicular sets of N parallel thin wires with Josephson 
junctions at each node (Figure 1) that is placed in an external tranverse field H. The classical thermodynamic variables 
of this array are the 27V superconducting phases associated with each wire. In the absence of an external field the phase 
differences would be zero at each junction, but this is not possible for finite H and the phases are thus frustrated. Here 
we assume that the Josephson couplings are sufficiently small so that the induced fields are negligible in comparison 
with H; we shall return to this topic when discussing the experimental realization of this network. Wc can therefore 
describe the array by the Hamiltonian 
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^ ^ ^'m^!Tmn'^n (1) 



where Jmn is the coupling matrix 



with Jjfe = exp{27riaj k / N) and 1 < (j, k)<N where j{k) is the index of the horizontal (vertical) wires; Sm = e*"^™ 
where the 0m are the superconducting phases of the 2N wires. Here we have introduced the flux per unit strip, 
a = NHP/(f)Q, where I is the inter- node spacing and 0o is the flux quantum; the normalization has been chosen so 
that To does not scale with N. 

Because every horizontal (vertical) wire is linked to every vertical (horizontal) wire, the number of nearest neighbors 
in this model is N and it is accessible to a mean-field treatment. For the same reason, the free energy barriers separating 
its low-temperature solutions scale with N. This situation is in marked contrast to that in conventional 2D arrays 
where the coordination number and hence the barriers are low. 1^,0 A similar long-range network with positional 
disorder was studied previously. For a 3> 1 /N the latter system displays a spin glass transition which was mapped 
onto the Sherrington-Kirkpatrick model Q for a 3> 1; in this field regime there is no residual "ferromagnetic" phase 
coherence between wires. Physically this glassy behavior occurs because the phase differences associated with the 
couplings Jjk acquire random values and fill the interval (0, 2tt) uniformly for a ^ 1/N . More speciflcally, there will 
be no commensurability if the sum 

EWL-Eexp27rz|^i^|fc (3) 

fc=l k=l ^ ^ 

is a smooth function of (j — /) where {j, V\ (k) are indices labelling horizontal (vertical) wires; this will occur only 
if the expression in curly brackets on the r.h.s. of (^ is not an integer for all k, a condition always satisfled for the 
disordered array. For the periodic case, this situation is realized in the "incommensurate window" 1/N <C a < 1; here 
the phase-ordering unit cell is larger than the system size so that the "crystalline" phase is inaccessible. ||^ There are 
thus no special fleld values where the number of low-temperature solutions are finite, in contrast to the situation for 
a > 1. 

In the thermodynamic limit of — s- cx) (with fixed array area), the high-temperature approach to the glass 
transition in the periodic array has been studied [^[^] using a modified Thouless-Anderson-Palmer (TAP) method. 

Here we discuss the qualitative picture that emerges from these results (cf. Fig. 1), referring the interested 
reader elsewhere for a more detailed quantitative treatment. As T approaches T+, where ~ Tq « there 
appear a number of metastable states in addition to the paramagnetic free-energy minimum. There is no dynamical 
transition at this temperature, indicating that the system is not "trapped" . We speculate that this is because these 
states are energetically unfavorable and/or they are separated by low barriers. At lower temperatures, as T ^ Tq , 
the paramagnetic minimum is "subdivided" into an extensive number of degenerate metastable states separated by 
effectively infinite barriers, and the system is dynamically localized into one of them. Qualitatively, in the interval 
Tm > T > Tq there appear many local minima in the vicinity of the paramagnetic state separated by Unite barriers; 
these barriers increase continuously and become infinite at T = Tq. Each of these minima is characterized by a 
finite "site magnetization" m.i = (sJt where"site" refers to a wire. When T > Tq thermal fluctuations average over 
many states so that {m.i) =0. Ai T — Tq the system is localized in one metastable state and there is an associated 
jump in the Edwards- Anderson order parameter, {q — jj "^ii^i)"^) ■ The low-temperature phase is characterized by 
a finite q and by the presence of a memory, limf^oo A(<, t') ^ where A{t, t') is the anomalous response. We expect 
that at T = Tq, the metastable states are degenerate and thus there can be no thermodynamic selection. However 
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at lower temperatures interactions will probably break this degeneracy and select a few states from this manifold; 
then we expect an {t oo) equilibrium first-order transition (Ts) which should be accompanied by a jump in the 
local magnetization. In order to observe this transition at Tg the array must be equilibrated on a time-scale {tw) 
longer than that {Ia) necessary to overcome the barriers separating its metastable states; tA scales exponentially with 
the number of wires in the array. Thus the equilibrium transition at Ts is observable only if t\Y — > oo before the 
thermodynamic limit (N oo) is taken; in the opposite order of limits only the dynamical transition occurs. 

The periodic array thus exhibits a first-order thermodynamic transition preceeded by a dynamical instability; the 
glass transition at Tq is characterized by a diverging relaxation time and an accompanying jump in the Edwards- 
Anderson order parameter. Furthermore the dynamical coupled equations for the response and correlation functions 
of this non-random network in the regime (1/iV <C a < 1) are identical to those obtained for the p = 4 (disordered) 
spherical Potts model. By contrast the behavior of the disordered array is similar to that of the p = 2 case 
with coinciding static and dynamic transitions. The possible connection between non-random glasses and p > 3 
(disordered) spherical models has been previously suggested in the literature, and the periodic array results 

give support to those conjectures. Furthermore these long-range Josephson arrays can be built in the laboratory, 
allowing for parallel theoretical and experimental studies of the "simplest spin glasses" that have previously been an 
abstraction. 

Since the experimental realization of these arrays is important, we would like to identify a number of simplifying 
assumptions inherent in the theoretical treatment discussed here; in order to test its predictions and probe beyond 
its realm (finite-size effects etc.) the fabricated system must satisfy certain physical requirements. In particular the 
model has been only considered in the classical limit where quantum mechanical fluctuations of the phase variables are 
negligible. Furthermore the fields induced by the Josephson currents must always be small compared to H so that it 
can frustrate the phases effectively. These constraints put some strong restrictions on the choice of array parameters, 
which we elaborate below. 

We begin with the first condition, that of maintaining weak quantum fluctuations. The strength of the latter is 
controlled by two parameters: the dimensionless normal resistance e^R/fi and the ratio Eq/Ej where Eq and Ej 
are the Coulomb and the Josephson energy scales respectively associated with each individual phase variables. The 
state of any array (or junction) is conveniently described by its position on the plane defined by these dimensionless 
parameters, Ec/Ej and e^R/fi; the resulting plot is called a Schmid diagram. ||T^ For a conventional short-range 
array R ^ Rq and Ec ~ e^/Co where Rq (Cq) is a resistance (capacitance) of an individual junction. We require that 
the capacitance of the individual wires and the dissipation in the system should be sufficiently large to suppress the 
quantum fluctuations of phase and charge effects. 

According to the Ambegaokar-Baratoff relation, the normal-state resistance of a single junction in the fabricated 
long-range arrays should be large (i?o S> ^/e'^) so that the associated screening current is small; for the sake of 
example we shall assume Rq w 100 kOhm and an individual capacitance of C « 1. 10~^^ F, parameters associated 
with preliminary Al — AI2O3 — Al tunnel-junction arrays already fabricated. j|] We note that the Josephson effect for 
a single tunnel-junction of such a large resistance (or for a short-range array of such junctions) would be completely 
overwhelmed by quantum fluctuations and Coulomb blockade effects. [pT|-po| Indeed, it is known that conventional 
nearest-neighbor arrays exhibit insulating behavior if their individual junction resistivitance is more than IbkOhm. 

However in the long-range Josephson arrays all the junctions associated with a particular wire are connected in 
parallel, so that the effective resistance and capacitance determining their position on the Schmid diagram ||l5[| are 
R — Rq/N and C — NCo respectively. Even for a small array with N — 100 the effective resistance and capacitance 
are R ^ 1 kOhm and C ^ 1. 10^ ^'^ F (corresponding to the charging energy Ec ^ 10 mK) respectively so that 
quantum fluctuations and charging effects are negligible for these networks. We note, however, that if the number 
of wires (and hence the number of neighbors) were decreased significantly the parameters of the individual junctions 
would have to be altered to avoid complications of charging effects. 

We would also like to minimize the effects of everpresent screening currents to ensure that the external field frustrates 
the superconducting phases effectively. More specifically we demand that the induced flux in the array is less than a 
flux quantum 



so that all fields and phase gradients produced by the diamagnetic currents are negligible. In order to determine the 
maximum induced flux in the array we need only look at its largest loop; our condition then becomes 



where I is the distance between nodes in the array. For the frustrated case, the maximum field occurs when the 
currents add incoherently in one direction and coherently in the other B so that 



Lieff < $0 



(4) 



(5) 
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ieff = (V^A^) ic- 



(6) 



where is the critical current associated with an individual junction. Using 



we find that the necessary condition for weak induced fields becomes 

Thus there is a delicate balance to be obtained between the need for many neighbors (and high free-energy barriers) 
and the requirement (^) that the effective London penetration length be substantial larger than the system size; 
unfortunately early experimental attempts to realize these arrays were not able to satisfy these restrictions due to 



resolution limitations associated with the optical lithography methods used. [£1 22 1 Furthermore the predicted glassy 
behavior is only observable in the field range (/)o/{Nl)'^ < H < 4>o/l'^; where the minimum and maximum correspond 
to a single flux quantum in the whole array and in a single plaquette. Even with the use of the e-beam lithography 
{I = 0.5 /im) the requirement (||) can only be satisfied for an array of Al — AI2O3 — Al junctions {Ej = 1. 10~^^ Erg) 
with the number of wires TV less than 300. For these parameters with N = 100 (j)o/{Nl)'^ is of the order of 0.01 Oe, 
and a careful screening of the earth's magnetic field is necessary. 

Once the arrays are fabricated to specification what measurements should be performed? The field-dependence of 
the critical temperature of both disordered and periodic arrays is crucial in establishing the effectiveness of of the 
external field in frustrating the phases. It is also a good test of "phase freezing" into a glassy phase, and the functional 
dependence Tc{H) cx is expected over a wide field and temperature range for both types of arrays. In order to 

firmly establish the presence of large barriers imperative for glassiness, the history-dependence of the critical current 
jc should be studied (i.e. its dependence on the path in the T — H plane) in both arrays. We note that below 
Tc(a) we expect a faster increase in jc(T) for the periodic network than for its disordered counterpart; |^ specifics 
for the former case have yet to be determined. We expect that the diverging relaxation time at Tc can be accessed 
experimentally via the a.c. response to a time- varying magnetic field H(t); the associated ac susceptibility is 

^ dM^ ^ C{a,N)uj 

dH^ uj + i/tR{T) ^ ' 

where tjf_ is the longest response of the system that diverges at T ~ Tq and C{a,N) is to be found elsewhere. ||] 
The cj — > limit of the a.c. susceptibility jumps to a finite value at T = Tq, indicating the development of a finite 
superconducting stiffness at the transition. Therefore measurement of this a.c. response in a fabricated array would 
probe its predicted glassiness. 

Work is currently in progress to study the physical properties of the periodic model in its non-ergodic regime 
(T < To). We note that the glass transition temperature Tq described above corresponds to the system going 
out of equilibrium as it is cooled infinitesmally slowly from high temperatures. In practice there will always be a 
finite-rate of cooling and the effective glass transition will occur when the system drops out of equilibrium. If we define 

the reduced temperature Q = (^^^^^~^^ then the effective glass transition will occur when the time associated with 
cooling, tc ~ dB/dt ^® equal to the relaxation time at that temperature, tji « tf)9^'' where to and are determined 
from the high-temperature analysis elsewhere; Q the effective glass transition temperature and the associated time are 
indicated as Tq and to in Figure 2. There is evidence for history-dependence; the system's response is very different 
if an additional field is turned on at a time, tn, during or after the cooling process and a subsequent measurement 
is taken at time tM (cf. Figure 2), reminiscent of the zero-field cooled vs. field-cooled susceptibility observed in spin 
glasses. Furthermore an "ageing" time-scale t* seems to exist in this system (cf. Fig. 2); if the field is turned on 
during the cooling process at time tn and a subsequent measurement is taken at time t^i such that tn — ^a/ < t* 
the system "remembers" the cooling process; otherwise {tn — tM > t*) it "forgets" it completely. Preliminary results 
suggest that t* has a decreasing temperature-dependence and probably vanishes at low temperatures. This presence 
of an additional "ageing" time-scale does not exist in the Sherrington-Kirkpatrick theory, though it is indeed observed 
in experimental spin glasses; it is thus amusing that it emerges from the study of a non-random model. 
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Fig. 2 Schematics of (a) Temperature (T) (b) Application of an Additional Field (Aif)vs time (t) indicating the time-scales involved 

in a finite-cooling experiment as discussed in the text; Tq and ta refer to the "effective glass transition" when the system goes out of 
equilibrium, and tu , tM and t* are the time-scales associated with the onset of an additional field, a subseqent measurement and ageing. 



In conclusion wc have presented a summary of recent work on a periodic model that displays "freezing" into one 
of an extensive number of metastable states without thermodynamic selection. This glass transition is characterized 
by a diverging relaxation time and an accompanying jump in the Edwards- Anderson order parameter. At low field 
strengths the dynamical behavior of this system is identical to that of the p = 4 (disordered) spherical model. A 
key advantage of this system is that it can be built in the laboratory, thereby allowing for the possibility of detailed 
interplay between theory and experiment. The behavior of a similar array with a reduced number of neighbors is 
a completely open question that could be important in understanding the relevance of long-range models to real 
experimental glasses. Finally, since any uncertainty in the position of the wires introduces randomness in the array, 
the continued study of this network could also offer an opportunity to study the crossover between glasses with 
spontaneously generated and quenched disorder. 
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